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ABSTRACT 

Building on a recent work of C. Crnkovic, M. Douglas and G. Moore, a study of multi- 
critical multi-cut one-matrix models and their associated s/(2,C) integrable hierarchies, 
is further pursued. The double scaling limits of hermitian matrix models with different 
scaling ansatze, lead, to the KdV hierarchy, to the modified KdV hierarchy and part of the 
non-linear Schrodinger hierarchy. Instead, the anti-hermitian matrix model, in the two-arc 
sector, results in the Zakharov-Shabat hierarchy, which contains both KdV and mKdV as 
reductions. For all the hierarchies, it is found that the Virasoro constraints act on the 
associated tau-functions. Whereas it is known that the ZS and KdV models lead to the 
Virasoro constraints of an s/(2,C) vacuum, we find that the mKdV model leads to the 
Virasoro constraints of a highest weight state with arbitrary conformal dimension. 
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Introduction 



Most of the interesting properties of the matrix model formulation of two dimensional 
gravity were originally extracted for the special case of lagrangians with even potentials [1] . 
In the hermitian 1-matrix model it was later shown that the introduction of odd terms in 
the potential gives rise to a doubling of the critical degrees of freedom and a doubling of the 
critical equations. In the 1-arc sector one gets two decoupled Painleve I equations, for the 
first critical point; the underlying integrable structure being two decoupled Korteweg-de 
Vries (KdV) hierarchies [2,3]. The situation, however, turns out to be different in the 2-arc 
sector of the theory, which, for an even potential, has a Painleve II equation as the lowest 
multi-critical point; the underlying integrable structure being the modified Korteweg-de 
Vries (mKdV) hierarchy [4,5]. Explicit calculations show that the introduction of odd 
terms do not lead to decoupled equations and to the doubling of the mKdV system [6] . 

This paper originated from our attempt to explore the integrable structures associated 
with the 2-arc sector of the hermitian 1-matrix model, with a generic potential. Indeed, one 
of the most interesting features of matrix models is the fact that the known 2d quantum 
gravity models (both pure and coupled to minimal conformal matter) are described by 
an integrable hierarchy, supplemented with an additional condition known as the 'string 
equation' [7]. The common belief is that the (anti-) hermitian n-matrix model should 
correspond to a hierarchy associated to the Lie algebra sl(n + 1,C), in the sense that, 
choosing different scaling ansatze for the double scaling limit of the matrix model, one 
gets a field theory, where the free energy (or a function related to it), satisfies the string 
equation of such a hierarchy. Obviously, the case n = 1 is the simplest and also virtually 
the only one where examples can be worked out explicitly. More general one-matrix models 
have also been considered in ref. [8], which deals with the case for complex matrices. 

We start by deriving in an explicit way the higher multi-critical points of the her- 
mitian 1-matrix model in the 2-arc sector with generic polynomials. We find that the 
resulting string equations are associated with only the 'even' subset of flows of the non- 
linear Schrodinger (NLS) hierarchy. This, indeed, seems to be the hierarchy behind the 
2-arc sector of the hermitian 1-matrix model, as one can check by computing explicitly the 
Lax operator. 

Interestingly, we notice that if the odd terms in the potential were purely imaginary, 
i.e. if the starting matrix model were anti-hermitian instead of hermitian, the multi- 
critical points in the 2-arc sector reproduce all the string equations associated to the ZS 
hierarchy. Indeed the Lax operator one gets in this case is that of the Zakharov-Shabat 
(ZS) hierarchy 1 . The reason why we do not find the odd string equations of the NLS in our 
derivation of multi-critical points of the hermitian matrix model is that those equations 
are complex. Instead all the flows of the ZS hierarchy are real as we will see in section 2.5. 

While the NLS hierarchy is known to contain, as reduction, the mKdV hierarchy, the 



The connexion of the ZS hierarchy with the (anti-) hermitian matrix model in the 2-arc sector 
was first noticed by the authors of ref. [9] . 
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ZS hierarchy contains both the mKdV and KdV hierarchies, [10], (which can equivalently 
be obtained directly from the matrix model by restricting the ansatz made for the double 
scaling limit). In addition, we find that the string equations of the NLS hierarchy reduce 
to those of the mKdV hierarchy, and the string equations of ZS to those of the KdV and 
mKdV hierarchies. From this one deduces that solutions of the KdV and mKdV theories 
correspond also to solutions of the NLS and ZS theories, in a particular subspace spanned 
by the even flows of the two hierarchies. The results are conveniently illustrated in diagram 
1, which shows how the critical points of the various hierarchies relate both in the 1-arc 
and 2-arc sectors. 

KdV 




even p. 



1-arc 

general p. 

KdV 

(double) 



NLS — ► mKdV 

(even) red. 



general p. 



.. even p. 

nerm. 



2-arc mKdV 

anti-h. 



even p. 

^ general p. KdV 
/ 

ZS red. 

\ 

mKdV 

Diagram 1 

We then go on to discuss the tau-function formalism of the ZS hierarchy. In this 
case we find that the partition function of the theory is equal to the tau-function of the 
hierarchy, rather than to its square, as happens instead in the 1-arc KdV model. As it 
was first shown in ref. [9], the partition function of the matrix model leading to the ZS 
hierarchy, satisfies the Virasoro constraints of an untwisted boson with an, a priori, arbi- 
trary value for the zero-mode or 'momentum'. This implies that the Virasoro constraints 
act on the tau-function, to mirror the situation for the 1-arc KdV model [11]. However, 
in the ZS case, the tau-function carries an additional quantum number due to the zero- 
mode of the untwisted field. This additional quantum number seems to play the role of 
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a non-perturbative parameter which labels different sectors of the theory and arguments 
connected with the tau-function formalism suggest that it takes discrete values. 

An interesting side issue concerns the existence of Virasoro constraints for the models 
described by the mKdV hierarchy; these include both the 2-arc (anti-) hermitian model, 
with even potential, and the unitary matrix models [12,13]. We find that there are, indeed, 
Virasoro constraints for the mKdV model; however, the situation is more complicated than 
in the KdV case. In the KdV case, the Virasoro constraints act on the tau-function: that is 
the square root of the partition function. For the mKdV model there are two tau-functions; 
the partition function being the product. We find that the string equation and the mKdV 
hierarchy imply a set of Virasoro constraints for each tau-function separately, however, in 
contrast to the KdV case, only the L n for n > appear, and the eigenvalue of L is an, 
a priori, undetermined integration constant. In other words, the tau-functions satisfy the 
Virasoro constraints of a highest weight state of the conformal algebra. 

Under the reduction from ZS to KdV, the Virasoro constraints are transformed into 
those of a twisted boson acting on the square root of the partition function (as expected 
[11]). More interestingly, the reduction to KdV of the ti = constant, to = x, ZS scaling 
theory gives rise to the topological point of the KdV hierarchy! In other words, from the 
anti-hermitian 1-matrix model with a potential of the fourth order it is possible, after the 
double scaling limit, to get Topological Gravity [14]. 

For the reduction to mKdV we get a series of constraints acting on the mKdV partition 
function which are not Virasoro-like; this is to be expected because they act on the product 
of the mKdV tau-functions and not on each of them separately. We also show that the 
'conformal dimension' of the mKdV partition function is fixed under the reduction. Finally 
in both reductions the value of the above-mentioned discrete parameter is fixed. 

Another important issue which arises from this analysis, is the possible existence of 
new continuum theories described by the 'odd' critical points of the ZS model. The first 
such model was discussed in [9], and it seems to display a topological nature, similar to 
the topological point in the 1-arc KdV model [14]. Nevertheless, the theory is sufficiently 
different from the 'normal' topological theory, to make it unclear as to what continuum 
theory it describes. However, like its 1-arc cousin, this topological theory cannot be ob- 
tained from the matrix model. The possible higher 'odd' multi-critical points have yet to 
be explored. We will not address the mathematical technicalities required to prove that 
solutions of these new scaling points exist, however, since the mathematical apparatus 
exists [15] we hope these questions will be tackled elsewhere. 

Finally we notice that models with anti-hermitian matrices seem to arise in relation 
with topological gravity (the Kontsevich model [16]) and with the Penner model [17]. 
Indeed, whereas hermitian matrix models naturally emerge from the study of quantum 
gravity as a theory of random surfaces [1] , anti-hermitian matrices seem to arise when one 
tries to make connexions between matrix models and moduli spaces of Riemann surfaces. 



1. (Anti-) Hermitian 1-Matrix Models and Orthogonal Poly- 
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nomials 



In this section we will consider the double scaling limit of a general (anti-) hermitian 1- 
matrix model. Most of what follows has been already extensively discussed in the literature, 
so we will recall only the most important results, add some new ones and develop a few 
explicit examples to fix our notation. 



1.1 Double Scaling Limit and String Equations of Hermitian Matrix 
Models 



Let M be a hermitian N x N matrix and consider 



'N 



= j dMe-^ v{M ^ (1.1) 



where V(X) = g±X + ^-A 2 + A 3 + • • • and A denotes a (real) eigenvalue of M. 

As it is well known, one can introduce orthonormal polynomials V n (X) such that 



J dXe-? v ^V n (X)V m (X) = 5 n , m (1.2) 



with 



XV n (X) = v ^P n+1 (A) + 5 n ? (l (A) + y/RnV n -l(X) , (1.3) 

then 

N-l N-l 

Z N = M H hi = M C n R?-* (1.4) 

i=0 i=l 

where R n = h n /h n -\. 

We will consider, for the moment, the case with real potentials. Notice that this 
implies that both R n and S n are real. Moreover, if V(X) = V{— X) then S n = 0. 

The most important equation is the so-called 'string equation', which can be written 

as: 

Tl 

"5 = <72-Rn + gz{RnS n + RnSn-l) H 

P (1.5) 
= 9i+92S n + g 3 (Sl + R n + R n+1 ) + -- - . 

In the double scaling limit one assumes that j3/N — > 1, n/(3 = 1 — x/N a and 

_ , ^s n f(x) r(x) j(x) z(x) 
R n = ! + {-!) — + — + — + — + ■■■ 

C - hM( -I\n9(x) , S(X) p(x) V(X) 

where b is an arbitrary constant. 
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For an even potential, V(X) = V(-X), and V(X) = fX 2 H + A 2fc , one sets in 

the 1-arc sector 

2k i 
f(x) = = g(x), «=2fcTl' 7i = ^ = 2fcTl (L7) 

and in the 2-arc sector 

2fc~2 % 

a= 2i—T> 74 = * = 2fc=T ' (L8) 

and one needs to introduce scaling functions up to the order 2k — 2. 
Let us introduce a such that 2 

a = N^k+r in the 1-arc sector 
i 

a = A^ 2fe -! in the 2-arc sector, 

then eqs. (1.5) become 

x f { Xl i 
1 = Fo(9i) H ^i(^) + — F2(f,9,r, s, f',g';gi) + ■■■ 

= ^o(^) + ^^i(^) + ^2(/^,r, S ,/ / ,^;^) + --- 
a ar 



(1.9) 



(1.10) 



where e = 2/c in the 1-arc sector, e = 2/c — 2 in the 2-arc sector and ' means 4-. The 
functions JF^ and ^ can be explicitly obtained from eqs. (1.5) and depend on the coupling 
constants and on the scaling functions. 

It is now easy to solve order by order in 1/a the string equation (1.10), for example 
the zeroth order in 1/a (JF = = Go) fixes the value of g\ and g^\ 

o 2 = 1 - 26a 3 - (3 + 3b 2 )g 4 - (12b + 4b 3 )g 5 - (10 + 306 2 + 5b 4 )g G + ■■■ 
9i = - bg 2 - (2 + b 2 )g 3 - (66 + b 3 )g 4 - (6 + 126 2 + b 4 )g 5 + (1.11) 
- (306 + 206 3 + 6 5 )a 6 + --- • 

Notice that the g^s depend on the free parameter b. 

At order 1/a one finds two possible solutions: either f(x) = and g(x) — or one 
fixes the value of g 4 . Choosing the solution f(x) = = g(x) one gets the 1-arc sector 
string equations; fixing g 4 instead, one gets the 2-arc string equations. In the case of an 
even potential of order 4 (k = 2) these are the well-known Painleve I and Painleve II string 
equations. 

Moreover, from eq. (1.4), in the double scaling limit, we find the 'specific heat' (up 
to non-universal terms) 

F" = \f 2 {x)-r{x) (1.12) 



2 This is not the 'a' defined by Na 2+1 ^ m = 1 usually introduced in the literature. 

3 Most of the computations of this section have been done with the help of Mathematical^ 1 . 
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where Z = exp(-F) and (PP) = <9 2 \ogZ = —F". 

Let us consider now a real potential V(X) with both even and odd couplings. In 
the 1-arc sector (f(x) = = g(x)) one gets the well known 'doubling' phenomenon. For 
example, from 

V(A) = (±& 3 -26)A + ±(2-6 2 )A 2 + ±&A 3 -^A 4 (1.13) 
one has the string equations 

x = \\" + X 2 , x = \x" + X 2 (1.14) 

where X(x) = r(x)—s(x) and X~(x) = r(x)+s(x). Notice that these equations do not depend 
on the parameter b and that the (721+1 are proportional to b. In the same way, in the 1-arc 
sector one obtains all the string equations of the KdV hierarchy, which are described in the 
following chapter. The 'specific heat' for these models is given by F" = —r(x) = — \{x+X)- 
The 2-arc sector with a general real potential is more interesting, indeed many different 
scaling solutions and string equations appear. For example, with 

V{\) = (26-6 3 )A+ ±(36 2 -2)A 2 -6A 3 + ±A 4 (1.15) 

(1.16) 



one gets 



f"-\f{g 2 + f 2 ) + \fx = 
9" - \g(g 2 + f) + \gx = o 

and from 

V(A) = (b + b 3 - |6 5 )A + |(-1 - 36 2 + |6 4 )A 2 + §(36 - 56 3 )A 3 + 
i(-l + 56 2 )A 4 -i6A 5 + ^A 6 

= xf - §/ V + f 2 ) - lfg 2 (f 2 + g 2 ) + lf(f') 2 + Sgf'g'+ 

- \f{g') 2 + Iff" + %g 2 f" + fgg" - / (4) 

=xg- §2 V + f) - lgf(f + g*) + lg(g'f + 3ff'g'+ 

- Uf? + \g 2 g" + lf 2 g" + fgf" - g [A) • 



(1.17) 



one gets 



(1.18) 



Again these equations do not depend on the parameter b; the (721+1 are proportional to 
b and putting g = these equations will turn out to be the first two string equations of 
the mKdV hierarchy [5]. When g is not put to zero, these two string equations will be 
seen, in the next chapter, to be the second and fourth string equation of the Non-Linear- 
Schrodinger (NLS) hierarchy. It is not possible to get the first and third string equation of 
the NLS hierarchy from the hermitian 1-matrix model. It seems that the NLS hierarchy 
does not have a full 'matrix model' realization, in the sense that we do not get the multi- 
critical points corresponding to the 'odd' flows; a fact that will be explained in section 
§2.5. 

For all these models the 'specific heat' turns out to be 



F" = \ [f\x)+g\x)] (1.19) 



since r(x) = j(f 2 (x) - g 2 (x)). 
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1.2 



Anti-Hermitian Matrix Models and Their String Equations 



Let us first consider, in more detail, the NLS string equations (1.16) and (1.18). 
Sending g — > ig one obtains string equations that are nothing but the string equations 
associated to the Zakharov-Shabat (ZS) hierarchy, as will be apparent from the next chap- 
ter. But this substitution is incompatible with our ansatz, eq. (1.6), for the double scaling 
limit. Indeed, sending g — > ig implies that S n becomes complex. But if the potential 
is real, V*(X) = V(A), it is easy to show that S n must be real. On the other hand, if 
V*(X) = V(—X) then S n is pure imaginary. (This can be easily shown using the orthogo- 
nality of the polynomials P n (X) (P n (X) = VKV n (X)), the fact that P n (A) = A n + 0(A n_1 ) 
and the reality of the partition function (see (1.4)).) Therefore, we are led to consider the 
potential V(A) = ig\X + \g2X 2 + |i#3A 3 + ^g^X 4 + ■ ■ ■ with A, e R and some of the 
g2i+i different from zero. The double scaling limit is now done in exactly the same way 
as in the previous section except for the ansatz for S n (eq. (1.6)) which is instead pure 
imaginary, i.e. 

s n = i b+{ -i r m + ^ + ^ + ^ + ... (1 . 20) 

a a 2 a 6 a 4 

with 6, g(x), s(x), ... all real. 

Notice that we can get rid of the l V in the potential introducing anti-hermitian ma- 
trices M = %M with pure imaginary eigenvalues A = iX. Again the R n are real and in the 
same way it is easy to prove that the S n are pure imaginary. 

For notational simplicity in the rest of this section we will continue to use hermitian 
matrices with complex potentials instead of using directly anti-hermitian matrices. 

In the 1-arc sector one gets again the double KdV string equation, now however in 
the variables X = r — is and X = r + is = X* ■ 

In the 2-arc sector we found the family of scaling solutions corresponding to the ZS 
hierarchy. For example, from 

V(X) = i(2b + b 3 )X + i(-2-36 2 )A 2 -ibX 3 + ±A 4 (1.21) 

with e = 2 one gets 



= \xf + f" + \f{g 2 -f) 

= l X g + g" + Igtf - f) . 

From 

V(X) =i[-l + b + b 2 -b 3 + \b 4 - i& 5 ] A+ 
i [_i_26 + 36 2 -26 3 + |6 4 ] A 2 + 
\i [-1 + 36 - 36 2 + 56 3 ] A 3 + \ [-1 + 2b - 5b 2 } A 4 + 
\i [J - §6] A 5 + ^ 



(1.22) 



(1.23) 



and e = 3 one has 



= xf + g"' + lg'(g 2 -f 2 ) 

= xg + r + lf(g 2 -f 2 ) . [ • } 
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Finally, from 

V(X) =i[b-b 3 - ±6 5 ] A + \ [-1 + 36 2 + |6 4 ] A 2 + \i [3b + 56 3 ] A 3 + 

i [_i _ 56 2 ] A 4 - \ ib\ 5 + ^A 6 " (L25) 



with e = 4 one has 



= x/ + |/ 3 (s 2 - f) - Ifg 2 (g 2 - f) + |/(f) 2 " 3(//V+ 

§/(</) 2 + I/ 2 /" - §<7 2 /" " to" - / (4) 
= xg - lg 3 (g 2 - f 2 ) + Igf 2 (g 2 - f 2 ) - lg(g') 2 + Sff'g'+ 

-Uf'?-b 2 9" + lf 2 9" + f9f"-9 {A) ■ 



These are the first three string equations of the ZS hierarchy (excluding the topological- 
like point [9] which cannot be obtained from the matrix model). More precisely, they 
correspond to the points t e = constant, to = x and U = for % ^ {0, e} in the ZS hierarchy 
(see next section). 

A few comments are in order. As usual the string equations do not depend on b. The 
<72«+i are proportional to b except for the case of eq. (1.23). Indeed, it is not possible to 
get the string equation eq. (1.24) from a real potential (V*(A) = V(A)). Moreover, setting 
(7 = one obtains the first two string equations of the mKdV hierarchy (eq. (1.24) vanishes 
identically) and setting \ (f + g) = —1 in eqs. (1.22) and (1.26) one gets the first two string 
equations (topological point included) of the KdV hierarchy in the variable ip = ^(f ~ d)- 

For all these models it turns out that 



F" = -\[g\x)-f 2 {x)] (1.27) 

since r(x) = \[f 2 {x) + g 2 {x)]. 



1.3 Lax Operators from Matrix Models 

The basic idea, due to Douglas [7] , is that it is possible to make the double scaling limit 
not only on the string equations but also directly on A and §^ seen as operators acting on 
^(A). This operator is then related to the Lax operator of the integrable hierarchy which 
underlies the behaviour of the continuum theory associated to the particular sequence of 
scaling ansatze chosen. Indeed, under a double scaling limit from eq. (1.3) it is easy to see 
that A becomes a differential operator (A) of order 2 in x. The obvious relation 

= 1 (1.28) 

after the double scaling limit becomes the string equation [7]. Eq. (1.3) under the double 
scaling limit becomes 

A* = A* (1.29) 



d_ 



A 
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where which can be a vector, is related to the rescaled polynomials and A is a differential 
operator of degree 2 in x. Analogously, there exists an equation of the form 



— = . (1.30) 

OA 



These two equations can be rewritten as 

def 




(A-A)* = 0, — -mU = (1.31) 



and then the string equation becomes the compatibility condition for these differential 
equations [15], i.e. 

"d 



ax~ M 



. (1.32) 



It is easy to get the explicit expression of L from eq. (1.29). L is then the Lax operator 
of the corresponding hierarchy. Using the techniques of Zakharov-Shabat and Drinfeld- 
Sokolov [18,19], one can also explicitly compute M [7,4,9]. 

We now explicitly compute L from eq. (1.3). Consider first the case of a hermitian 
matrix with a real potential. Let n(x, A) denote V n (X) after the double scaling limit 4 , thus 
eq. (1.3) becomes 

XV n (X) ~ (2 + b)U(x) + \[r(x)U(x) + s(x)U(x) + U"(x)] + --- , (1.33) 

or 

where for notational simplicity we have dropped the dependence on A in IT. Setting A — 
(b + 2) — > A/a 2 one has 

LU(x,X) = , (1.34) 

where 

L = <9 2 + (r(x) + s(x)) - A . (1.35) 

This is the Lax operator of the KdV hierarchy. 

Notice that setting U(x, A) ~ (-l) n P n (A) one gets L = &l + (r(x) - s(x)) + A. Thus 
there are two KdV Lax operators associated to the hermitian 1-matrix model in the 1-arc 
sector, this is nothing but the doubling phenomenon [2]. 

In a similarly way, for the anti-hermitian 1-matrix model in the 1-arc sector one gets 
the previous two Lax operators with s(x) — > is(x). 

Consider now the case of the 2-arc sector with a hermitian matrix and a real potential. 
In order to get a non trivial Lax operator in the double scaling limit, following [4], we 
need to introduce two scaling functions depending on parity, I1(A, x) ~ (— l) m "P2m(A) and 
fl(X,x- I) ~ (-l) m V 2 m+i(X) where x = x 2m - Thus for n even eq. (1.3) becomes 

A(-l) n / 2 P n (A) ~ bU(x) + - [g(x)U(x) - f (x)fl(x) - 2Q' (x)} + ■ ■ ■ (1.36) 

ct 



4 Actually, as explained in refs. [7,4,15], under the double scaling limit Il(:r, A) 
e~ NV ( x ^ 2 V n (X) has a smooth behaviour and should be considered for this computation. 
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and for n odd one has 

A (-l) (n -^/ 2 P n (A) ~ M(x) + - [-g(x)n(x) - f(x)U(x) + 2U'(x)} + ■■■ . (1.37) 



a 



Thus, setting A — b — > A/a one has 
where 



L^> = (1.38) 



^-^(i- i)^(!i) + s(- ii) (139) 

and $ = Now, after having conjugated L and rotated we get eq. (1.38) with 

L = ia 3 \ + a 2 g -cnf + d x (1.40) 

and [07, (Tj] = ieijk<?k- As we expected, this is the Lax operator of the NLS hierarchy. 

In the case of an anti-hermitian matrix with pure imaginary eigenvalues A one can 
do the same computation with V2m+i and S n pure imaginary (see eq. (1.20)). Letting 
A — ib — > iX/a, one finally gets 

L = d x + ia 2 g - aif - a 3 X , (1-41) 

and this is the Lax operator of the ZS hierarchy. 

Notice that the mKdV string equations can be obtained both from the NLS (hermitian 
matrix, real potential) and ZS (anti-hermitian matrix, complex potential) hierarchy. This 
is obvious since setting S n = (b = = g = . . .) hermitian and anti-hermitian matrix 
models coincide. 



2. Hierarchies and String Equations 

In this section we will discuss the hierarchies of integrable equations which lie behind 
the non-perturbative structure of the one matrix model. The relevant hierarchies are well- 
known in the mathematical physics literature, and are the Zakharov-Shabat (ZS), for the 
Lie algebra sl(2, C), the Korteweg-de Vries (KdV), modified Korteweg-de Vries (mKdV) 
and non-linear Schrodinger (NLS) hierarchies [10]. 

Interestingly, all these hierarchies are intimately related to the Lie algebra sl(2, C). 
This is apparent in the matrix Lax formalism [18] and also the 'Hirota' or 'tau-function' 
formalism [20]. It is also significant that the complex-ZS hierarchy is the 'master' hierar- 
chy for sl(2, C), since the ZS, NLS, KdV and mKdV hierarchies can all be obtained by 
appropriate reductions. This will be explained in more detail below. We shall also discuss 
how the hierarchies are related to the matrix models. 
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2.1 



The sl(2, C) Hierarchies 



In section §1.3 we showed how the matrix models are connected with the integrable 
hierarchies presented in the Lax formalism. In fact, the most economical way of explicitly 
introducing the hierarchies is via their recursion relations. These can be extracted from 
the Lax formalism, see for example ref. [18]. In what follows we often write x for the 
flow t , °f the hierarchy under discussion. For the KdV equation the hierarchy can be 
presented, in the following way 

Ou 

— = d x R k+1 k>0, (2.1) 
ot k 

where R k is a polynomial in u and its ^-derivatives (the R k s are known as the Gel'fand- 
Diki polynomials [21], and are not to be confused with the R k used in the previous section). 
The integrability of the hierarchy is summed up in the recursion relation 

d x R k +i = {\dl + 2ud x + u') R k . (2.2) 

Specifying Rq = 1 along with the recursion relation completely determines the hierar- 
chy. To make contact with the notation used in the previous section, for example for the 
hermitian matrix model in the 1-arc sector with even potential, one should set u = r. 

For the complex-ZS hierarchy a similar structure of recursion relations is found. For 
two independent complex variables ip and ip, if 

= \{F k+1 - Gk+i) , t^- = \{F k +i + Gk+i) , (2.3) 

with k > — 1, then 

F k+1 = G' k + (^-ij)H k 

G k+1 = F' k + $ + ^)H k (2.4) 
H' k =^(G k -F k )-i;(G k +F k ) . 

Specifying Fo = ip — tfj and Gq = ip + ip along with the recursion relations then determines 
the whole hierarchy. Notice that the flow t-i is particularly simple 



implying the following dependence on £_i: ip ~ e _t_1 and if) ~ e* _1 . To make contact with 
the notation of the previous section ip = \{f — g) and ip = |(/ + g), which also agrees 
with the conventions of ref. [9] when t k — > t k +\. The NLS hierarchy is recovered from 
the complex-ZS hierarchy by choosing tp to be the complex conjugate of ip (up to a term 
involving t_i), a choice which is easily seen to be consistent with the recursion relations 
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of the hierarchy. The ZS hierarchy itself, simply corresponds to the complex-ZS hierarchy 
with t(j and ip both real. 

The mKdV hierarchy is best introduced through its relation to the KdV hierarchy via 
the Miura Map. This map takes a solution of the mKdV v hierarchy into a solution of the 
KdV u hierarchy as 

u = -v' - v 2 . (2.6) 

The pull-back of the KdV flow tf dV under the Miura Map is then the mKdV flow tf KdY . 
We now make this more explicit. Defining D = —d x — 1v and D* = d x — 2u, the recursion 
relation (2.2) becomes d x Rk+i = D( — ^d x )D*Rk [22]. Since d x u = Dd x v, we have 

g^W = D g^W = D H*>*) > ( 2 ' 7 ) 

and so the mKdV flows are 

^ = -±d x D*R k k>0, (2.8) 
ot k 

where the polynomial R k = Rk\—v' — v 2 ] is expressed in terms of v and its x-derivatives 
by substituting u for v via the Miura Map. To make contact with the notation of the 
previous section v = f /2. 

2.2 The ZS, NLS, KDV and MKDV Hierarchies as Reductions of the 
complex-ZS Hierarchy 

In this section we explain how the the various hierarchies that we have introduced can 
all be obtained as reductions of the complex-ZS hierarchy. The situation is conveniently 
summarized by diagram 2 which shows how the various hierarchies are related by reduction. 



NLS — ► mKdV 



complex-ZS 

\j KdV 

/ 

ZS 

\ 

mKdV 

Diagram 2 

We have already noted how the complex-ZS hierarchy is reduced to the ZS and NLS 
hierarchies by taking two different 'real slices'. For the former one takes ip and to be 
real, whilst for the latter one takes the complex conjugate of ip to be ip* = e _2 *- 1 V'- 
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The KdV hierarchy is obtained from the ZS hierarchy by setting 



^ = -e l ~\ (2.9) 

the KdV variable being given by u = —iftip. One can readily prove that F 2k +i + G 2k +i = 0, 
when evaluated at ip = — e* _1 , whereas F 2k + G 2k 7^ 0. This means that only the even 
flows preserve the condition (2.9), and so only the even flows reduce to flows of the KdV 
hierarchy. One finds 

d(~^) 1 TTl I (t) im 



<%2fc 

where * means 'evaluate at ip = — e* -1 '. If one now compares the recursion relations of the 
ZS hierarchy (2.4), evaluated at ip = — e t_1 , to those of the KdV hierarchy (2.2), then one 
deduces that 

H 2k+ i[^ = -e*- 1 ] = 2 k+1 R k+1 [-^}. (2.11) 

This means that the relation between the flow variables is t 2k = 2~ k tf dV . It is straight- 
forward to see that the KdV hierarchy cannot be obtained from the NLS hierarchy by a 
similar reduction. 

The mKdV hierarchy is obtained from both the NLS and ZS hierarchies by setting 

^=e 2t - 1 il;, (2.12) 

the mKdV variable being given by u 2 = ipip. One finds that G 2k+ i = F 2k = F[ 2k = 0, 
when evaluated at (2.12). So the situation is similar to that for the KdV reduction, in 
that only the even flows preserve the reduction. By pursuing a similar analysis of the 
recursion relations, one discovers that the flow variables are related via tf fc s = 2~ k tf KdV 
since F 2k+1 = —2 k d x D*R k {—i>' — z/ 2 ], when evaluated at (2.12). 



2.3 The String Equations and Their Reductions 

The hierarchies that we have discussed above admit many types of solution. However, 
in applications to matrix models and two-dimensional field theories, very particular solu- 
tions are required. In addition to boundary conditions, these are specified by adjoining to 
the hierarchy an extra condition called the string equation [7]. The string equation must 
be consistent with the flows of the hierarchy, in the sense that it must be preserved by 
the flows of the hierarchy. It turns out that the string equation admits certain scaling, or 
multi-critical solutions. It is these solutions which are found in the matrix model, after 
the double scaling limit. They are obtained by restricting to the subspace t k = constant, 
to = x and tj = otherwise, for some k. We call the resulting reduced equation the k th 
string equation. The string equations can be found in general following the analysis of 
§1.3. 

The string equation associated to the KdV hierarchy was found originally in [7]. In 
our conventions it takes the form 

J2(2k + l)t k ^- = -1 , (2.13) 
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which may be integrated using (2.1) to give 



J2(2k + l)t k R k = . (2.14) 

fc=0 

The k th multi-critical point corresponding to t k = constant, to = and tj = otherwise, 
is described by the string equation 

(2k + l)t k R k [u] = -x . (2.15) 

For the ZS hierarchy, the string equation was found in [9] 

00 P> i 00 pTT 

k=o k 1 k=o k 1 

The k th multi-critical point for which t k =constant, to = x and tj = otherwise, is 
described by the string equation 

(k + l)t k (F k -G k ) =2xi;, (k + l)t k (F k + G k ) =2xi) . (2.17) 

The above equations also apply to the NLS hierarchy by taking tfj* = e~ 

The string equation of the mKdV hierarchy has been obtained in [4], however, one 
can obtain it in a simple way given the string equation of the KdV hierarchy by using the 
Miura map. The idea is to pull back the string equation of the KdV hierarchy via the 
Miura map u = — v' — u 2 ; the result is then guaranteed to be consistent with the flows of 
the mKdV hierarchy, because of the Hamiltonian property of the Miura map. We act on 
(2.14) with the operator ^d x D* and use eq. (2.8) to obtain 

OO 

^(2A; + l)t fc -f = -xv' - v . (2.18) 



dt k 
fc=i 



This can be rewritten as 



00 P> 

J2(2k + l)t k ^ + is = . (2.19) 

fc=0 k 

Using eq. (2.8) again, we can integrate with respect to x (discarding, an integration 
constant) obtaining the mKdV string equation 

J2( 2k + 1 )tkD*R k = 0. (2.20) 

k=0 

where R k = R k [— v' — u 2 ]. The k th multi-critical point for which t k = constant, to = x 
and otherwise tj = 0, is described by the string equation 

(2k + l)t k D*R k [-!/' - v 2 \ = 2xu . (2.21) 
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We now show that the string equation of the ZS hierarchy consistently reduces to the 
string equations of the KdV and mKdV hierarchies, respectively, for the reductions we 
discussed in §2.2. For the reduction to KdV we evaluate (2.16) at ip = — e*" 1 , t2k+i = 
and use the relations obtained in §2.2 to get 



^+I)^ dV 2- fc (^ + ^)^-i = 

k Z° ( 2 - 22 ) 



k=0 



These equations are obviously equivalent to the KdV string equations (2.13) since H 2 k-i = 
2 k R k at^= -e*- 1 , t 2k +i =0. 

Analogously, in the case of the mKdV reduction with ip = e 2 * -1 ^, £2/0+1 = 0, one of 
the equations (2.16) is trivially solved and the other equation becomes 

00 

5> + ±)C KdV 2" fc G 2fc = 0. (2.23) 

k=0 

But G-ik = -2 k D* R k [-u' - v 2 ] at ip = e 2 *- 1 -^, t 2k+1 = 0, which gives the mKdV string 
equation (2.20). 

The significance of the fact that the string equation of the ZS model reduces to that 
of the KdV and mKdV models, is that solutions of the KdV and mKdV string equations 
are, by pulling back, solutions of the ZS string equation (on the subspace t 2k +i = V/c). 
This will lead us to conclude that the model described by the ZS hierarchy includes the 
models described by the KdV and mKdV hierarchies. We further discuss these facts in 
§3.4 and §3.5. 



2.4 The Tau-Function Formalism 

There is an alternative formalism for constructing integrable hierarchies which was 
originally developed as a direct solution technique of the non-linear equations of the hi- 
erarchy. The central objects of this approach are the tau-functions, which satisfy a hier- 
archy of non-linear 'Hirota bilinear equations', see [20] for example. For the hierarchies 
that we are considering, the Hirota hierarchies are intimately related to the Lie algebra 
s/(2,C) (= Ai). In fact, they use the two vertex operator representations of the basic 
representations of the affine algebra A^ . This works in the following way. The basic 

representations of are carried by the Fock space of a scalar field, either twisted or 
untwisted. The Hirota equations for the tau-function are equivalent to the condition that 
the tau-function lies in the orbit of the highest weight state of the group associated to the 
affine algebra. The untwisted construction underlies the ZS and NLS hierarchy, whilst the 
twisted construction underlies the KdV and mKdV hierarchies [20] . 
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For the joint KdV and mKdV system (related by the Miura Map) there are two tau- 
functions To and t\ arising from the two basic representations of the Kac-Moody algebra 
The relationships between these and the functions u and v are 

u = 2d 2 x log to , v = d x log . (2.24) 

One of the equations of the hierarchy is 

rUn - 2t' t[ + t t[' = , (2.25) 

from which one can extract the Miura Map u = —v' — v 2 and the relation 

is 2 = -d 2 x \og (ran) . (2.26) 

The Hirota hierarchy which leads to the complex-ZS hierarchies has an infinite set 
of tau-functions. This is because the relevant vertex operator construction, in this case, 
involves an untwisted scalar field, which has a zero-mode. In order that the operator 
product expansions of the vertex operators are local, the zero-mode must be quantized, 
taking values in the weight lattice of the finite Lie algebra; A\ in this case, whose weight 
lattice is simply isomorphic to Z. We will label elements of the weight lattice with half- 
integers, so that the sub-lattice generated by the root consists of the integers. With this 
labelling, the integers and half-integers correspond to the two distinct basic representations 
of the affine algebra A^ . To make a connexion with the complex-ZS hierarchy one chooses 
a fixed element of the weight lattice, that is a half-integer n. Then 

^=^, ?=^, (2.27) 

where we have used a tilde in order to avoid confusion with the mKdV/KdV tau-functions. 
In addition, the equations of the hierarchy imply 

W=-d 2 x \ogf n . (2.28) 

The hierarchies for different choices of n are isomorphic. 

2.5 From Hierarchies to Matrix Models 

Given a hierarchy and its string equation, to obtain the field theory describing a 
matrix model after a double-scaling limit, one must first identify the partition function of 
the field theory with some variable in the hierarchy. 

If x is the cosmological constant then the 'specific heat' is F" = —d 2 log Z. Clearly 
— F" has a well defined scaling dimension, from a hierarchical point of view. So, in prin- 
ciple, one can construct all the terms of the correct dimension from the hierarchy, and 
then see which can be integrated twice with respect to x. However, even this would not 
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determine the normalization of the 'specific heat'. In the absence of any additional physi- 
cal requirements, one has to appeal to the matrix model. It transpires that the partition 
function of each particular model is related to the tau-function of the hierarchy in a simple 

way: 

ZzS = T n , ZKdV = Tq, ZmKdV = TqT\ , (2.29) 

where the tau-functions for the hierarchies where introduced in section §2.4. 

With x interpreted as the cosmological constant, the 'specific heat' is nothing but the 
correlation function —(PP) = log Z, where the operator which couples to the cosmo- 
logical constant is conventionally denoted P and called the puncture operator. Finally, the 
insertion of an operator in a correlation function is given by the corresponding flow in the 
hierarchy, for example 

(O t PP) = ^-(PP) (2.30) 

where Oq = P and the right-hand side is computed using the equation of the hierarchy. 
Moreover, with the identifications (2.29) 

(PP) ZS = -W = (PP) NLS , (PP) KdV = u , (PP) mKdV = -v 2 • (2.31) 

Let us consider first the hierarchies which should correspond to the anti-hermitian 
matrix models and their correlation functions constructed using the prescription above. 
The anti-hermitian 1-matrix model in the 1-arc sector should correspond to the sum of 
two KdV hierarchies in the variables X = r — is, X = r + is = X* ■ Although these variables 
are complex (r and s are real functions of x), it turns out that the correlation functions 
and the string equations are real. Indeed, the 'specific heat' is given by the sum of the 
'specific heats', 

(PP) = \{X + X) = r (2.32) 
and the correlation functions are given by 

(O k PP) = ±d x (R k+1 (x) + R k +i(x)) ■ (2.33) 

Using the recursion relations of the KdV hierarchy it is easy to show that Rk(x) = Rk{x)* 
and hence the correlation functions are real. Moreover, the sum and the difference of the 
string equations (2.14) for X and X give exactly eqs. (1.5) for this model. 

The anti-hermitian 1-matrix model in the 2-arc sector corresponds to the ZS hierarchy. 
First of all, notice that for the ZS and NLS models one has 

(O k P) = \H k+1 (2.34) 

implying, for example, (PP) = \H\ = —ipip. For the ZS hierarchy both ip and ip are real, 
and thus the whole hierarchy and all the physical quantities are real. The string equations 
(2.16) correspond to the sum and difference of the two equations (1.5) where the second 
equation is multiplied by i. 

We now turn our attention to the hermitian 1-matrix model. In the 1-arc sector we 
have a double KdV hierarchy where everything is expressed in terms of real functions. In 
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the 2-arc sector, instead, we found only some of the string equations of the NLS hierarchy. 
Indeed, for the NLS hierarchy, one can easily show that Fzk, Hik and G^fc+i are purely 
imaginary, whereas i^/c+i; -f^2fe+i an d Gik are real. Thus, although the 'specific heat' is 
real, many correlation functions are complex or pure imaginary. This obviously forbids an 
interpretation of the full NLS hierarchy as a field theory obtained after a double scaling 
limit of a hermitian 1-matrix model in the 2-arc sector. Taking the sum and the difference 
of the string equations (2.16) and using the recursion relations of the hierarchy one gets 
[9] 

oo oo 

Y,(k + l)t k G k = 0, ]T(£; + l)t fc F fc = 0. (2.35) 

fc=0 k=0 

The string equations (1.5) should correspond to the multi-critical points to = t n = 
constant and t k = otherwise. Thus, if n is even the first equation is real and the second 
is pure imaginary and they correspond to eq. (1.5) where the second equation is multiplied 
by i. These are the string equations we found in the first chapter. For n odd instead, both 
equations are complex and so they could not arise from the direct study of the matrix 
model. 



3. Virasoro Constraints 



In [11] it was shown, for the models described by the KdV hierarchy, that the string 
equation, along with the hierarchy equations, could be reformulated as an infinite number 
of Virasoro-like constraints on the square-root of the partition function of the model. These 
constraints have a natural interpretation in terms of the Schwinger-Dyson equations for the 
loops of the matrix model. For the KdV model the square-root of the partition function is 
the tau-function of the hierarchy. The fact that the Virasoro constraints act on the tau- 
function of the hierarchy seems to be a universal feature of all the models, as will become 
apparent. First we briefly review the case for the KdV hierarchy. Using (2.29), (2.24) and 
the string equation (2.13), and integrating twice we deduce 

£(fc+I)f fc — + to = 0. (3.1) 

k=l ' / 

Following [11], we use the recursion relations of the hierarchy (2.2) and the relation between 
To and u in (2.24), which together imply 

d 2 x -^—\ogT =(ldl + 2ud x + u')d x —\ogro, (3.2) 

in order to express the t k derivative in terms of the t k +i derivative. With this relation, 
one finds that L k To = implies L k +iTo = 0, where the L k are written below. All of the 
constants of integration encountered are set to zero, on the grounds that they would other- 
wise introduce spurious scales into the theory, except for the constant in the Lq constraint 
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which is dimensionless and fixed by the requirement that the algebra of constraints closes. 
The end result is that the tau-function satisfies an infinite number of constraints of the 
form 

L n ro = 0, n > -1 , (3.3) 

where the L n 's are the Virasoro generators of a Z2-twisted scalar field 

°° d 
L-i= E( m +^)^^ + l £ o 

m=l ° tm - 1 



L = ^(m+§)t m -— + ^ (3.4) 

m=0 OTm 

°° 9 n d 2 

L n =J2( m + ^ tm ai + i S 7* — ^ • 

m _g Vf'm+n m—1 m— 1 n— m 

3.1 Virasoro Constraints for the MKDV Hierarchy 

Although the matrix model which leads to the mKdV hierarchy has been discussed 
in the literature [4], the analogue of the Virasoro constraints do not seem to be have 
been determined before (although ref. [12] does discuss Virasoro constraints before taking 
the double scaling limit). In this section we find the constraints using the mKdV string 
equation (2.20) and the recursion relations for the hierarchy. 

The mKdV string equation (2.20) is obtained from integrating (2.19). From (2.19) 
one easily deduces 

00 ft 2 

J2(2k + l)t k ^ + 2v 2 = 0. (3.5) 

k=Q k 

Recall that the partition function of the mKdV model is equal to the product of the tau- 
functions To and t±. We now express v in (2.19) in terms of the tau-functions To and t±, 
using (2.24), and v 2 in (3.5) using (2.26). The resulting two equations can be decoupled 
to arrive at 



d. 



d 

^ (2/c + 1) ^ logTj 



X 

' k=0 



= , j = 0,1 . (3.6) 



Integrating twice, and eliminating dimensionful integration constants, the two resulting 
equations may be written simply as 

L Tj = HjTj , j = 0,1 , (3.7) 

where L is identical to the Virasoro constraint of the KdV model, eq. (3.4), and hq and 
Hi are two, a priori undetermined, dimensionless integration constants. They are not, 
however, independent as we now show. By substituting the expression for v in terms of 
the ratio t\/tq, in the equation for the flows (2.8), we deduce 



d_ 

dt k 



log (^j = -\D*R k . (3.8) 
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This can now be substituted directly in (2.20) to yield 

The above, along with eq. (3.7) implies that ^0 = ^1 = \x. 

Before we discuss the possible meaning of the constant fj,, we first present a simple 
argument for determining the higher Virasoro constraints. We already know from the 
construction of the KdV constraints, that LfcTo = implies Lk+iTo = 0. It is also straight- 
forward to verify that if LqTq = (xto then Litq = 0, regardless of the value of 11. Therefore 
we deduce that r satisfies the infinite set of constraints 

L n To = ixtq <5 n>0 , n>0. (3.10) 

To find the constraints satisfied by ti, we notice that T\ satisfies exactly the same recursion 
relation as tq, that is (3.2), except that u = 2<9^1ogro is replaced with u = 2<9^1og7~i. 
Therefore, the same arguments that were applied to determine the constraints on tq will 
be applicable to Ti, hence we deduce 

L n Tj = /iTj 5 nt o , n > , (3.11) 

for j = and 1. So the mKdV partition function is the product of two factors which 
separately satisfy a set of Virasoro constraints, however, in contrast to the KdV case there 
is no L-i constraint and the Lq constraint includes an undetermined constant. Notice that 
the requirement that the constraints form a closed algebra under commutation, does not 
in any way constrain the value of the constant. It is important to realize that To and T\ are 
not independent, in fact they satisfy a whole hierarchy of equations for which (2.25) is but 
the first. So although, at first sight, the mKdV Virasoro constraints look less restrictive 
than the KdV constraints, one must bear in mind the additional equations which tie To 
and T\ together. 

The appearance of a parameter, which is not determined from the matrix model, 
seems, at first sight, to be surprising. However, it is not totally unexpected, indeed, such 
an occurrence is found at the first critical point of the mKdV model. At this point, the 
square root of the specific heat, v, satisfies the Painleve II equation [4] 

v" - 2u 3 + xv = . (3.12) 

This equation is known to admit a one-parameter family of solutions [23]. The actual 
solution which describes the matrix model, requires a scaling behaviour u ~ z*, as x — > 00. 
Ref. [24] discusses how, for one particular value of the parameter, such a physical solution 
does exists and is unique. It would be natural to suggest that the parameter is related to 
fi, the eigenvalue of the L constraint. Indeed, (3.12) admits the trivial solution v = 0, 
which corresponds to the situation when fj, = j^, for which the Virasoro constraints have 
the solution r = 1 (i.e. r being the vacuum of the twisted Fock space). Notice that, for 
the KdV model, the Virasoro constraints are those of an sl(2, C) vacuum, whereas, for the 
mKdV model, the Virasoro constraints are those of a highest weight state of conformal 
dimension fx. 
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(3.9) 



3.2 



Virasoro Constraints for the ZS Hierarchy 



The analogous constraints for the ZS hierarchy and string equation were found in [9]. 
Here, we briefly repeat their derivation which leads to Virasoro type constraints for an 
untwisted boson. 

The string equations are (see eq. (2.16)) 

J2(k + l)t k F k = , J](A; + l)t fc G fc = 0. (3.13) 

fc>0 fc>0 

Consider first the objects 

Ii = J> + l)t k H' k+j = J> + l)tk(sG k +j ~ fFk +J ) ■ (3-14) 

fc>0 fc>0 

They can be reduced, using the hierarchy, to sums involving only F k , G k and their deriva- 
tives, which are related to the string equations. 

From Jo, integrating twice over x and introducing an arbitrary integration constant 5 
a, we get 

J> + l)f fc (O fc _ 1 ) + ^ = (3.15) 
fc>i 

which will lead to the L_i constraint. 
From Ii we get 

J2(k + l)tk(O k }+P = , (3.16) 

fc>0 

where we have picked up a new dimensionless integration constant, (3. This leads to the 
Lq constraint. 

Using a similar procedure, from I 2 we get 

J2(k + l)t k (O k+1 ) + a(P) = (3.17) 

fc>0 

which leads to the L\ constraint. 

Finally, with a few more steps from I 3 we obtain 

J> + l)t k (O k+2 ) + (P) 2 + (PP) + a(O t ) = . (3.18) 

fc>0 

Whilst the previous equations involve only first order derivatives of the partition function, 
the equation coming from I3 has second order terms, which fix the function on which the 
Virasoro constraints act. In fact, we can write 

(P) 2 + (PP) = (F'f - F" = \d\Z (3.19) 

5 We will discard all integration constants which would have non-trivial dimension. 
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and, therefore, the Virasoro constraints act on the partition function. Since for this model 
the partition function is equal to the tau-function, we find that the Virasoro constraints 
act on the tau-function, mirroring the situation for the KdV model. 

By consistency, the commutator of two constraints should be a new constraint on the 
partition function. Therefore, using [L n , Li] = (n — l)L n+ i with n > 2, we get an infinite 
set of constraints acting on the partition function. These constraint are the Virasoro 
constraints 

L nT zs = , n > -1 (3.20) 



where 



fe>i 

k>0 

V^V, d d 



(3.21) 



d ^ d 2 d 



Ln = J2(k + l)t k — + J2 rl „ — + a- 



/•>() ' dtk+n f^ Q dtkdt n -k-2 dt n -i 

and (3 = a 2 /4 has been fixed by the relation [Li, L-{\ = 2(Lq + a 2 /A — (3) = 2Lq. 

3.3 Connexion with the Tau-Function Formalism 

It was noticed in ref. [9] that the Virasoro constraints (3.21) are those of an untwisted 
scalar field. In the convention for which 

ip(z) = q — ip log z + i a n , (3.22) 

L — ' n 

neZ^O 

and [a m , a n ] = nd n+rjlj o and [q,p] = i, we have for k > 

V2 d 1 

t k = -j^a-k-i, Wk = ^a k+1 . (3.23) 

The zero- mode p is related to the integration constant a via p = a/y/2. The conjugate 
variable to p does not appear in the Virasoro operators. 

The fact that the Virasoro constraints are those of an untwisted scalar field is very 
natural from the point of view of the tau-function approach, which we explained in section 
§2.4 based on ref. [20]. Recall that the zero mode of the scalar field of the construction 
of ref. [20] must have the quantized values m/y/2, for m G Z, in order that the vertex 
operators have local expansions. The tau-function can be projected onto eigenspaces of 
the zero-mode, these are precisely the r n which where introduced in §2.4, with \/2n being 
the eigenvalue of the zero- mode. 
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From eq. (2.29) the partition function of the ZS model is equal to f n , for some fixed 
half-integer n. It is very natural to identify the scalar field of the Virasoro constraints 
with the scalar field of the Hirota equations of [20]. We do not yet have a direct proof 
of this, however, below we present some arguments which support this view. Given this 
identification, one is led to a relation between the parameter n and the integration constant 
of the Virasoro constraints a: 

p = = -V2n, n e ±Z . (3.24) 
v2 

The possibility that a is quantized seems to be consistent with the results that we obtain 
in §3.4 and §3.5, for the KdV and mKdV reductions which require a = and a = — 1, 
respectively 6 . To substantiate this identification we now show that if f n satisfies the 
Virasoro constraints eq. (3.21) with a, then T n ±\ satisfy the same constraints but with 
a — > a=F2. We prove this fact following a similar demonstration as the previous paragraph. 
Let us consider the objects 

Y 3 = J2(k + l)t k -?—( f -f±-) . (3.25) 

and use the string equations, the hierarchy of Hirota equations for the f n 's, and the fact 
that f n satisfies Virasoro constraints with a, L m f n = L m (a)f n = 0, m > — 1. Considering 
Yq, Y\ and Y 2 , it is easy to obtain 

L_i(a =F 2)f„±i = L (a T 2)f n±1 = L x {a T 2)f„±i = . (3.26) 

Again, the L 2 constraint is more tricky. Considering Y" 3 , it is easy to show that 

d _ d 2 \ ( f n±1 \ 2 



L 2 (a T 2)f n±1 = (a T 2)|f n ^— J T 2f n±1 ^log f n j (3.27) 

The right hand side of this equation vanishes because of one of the Hirota equations satisfied 
by the tau-functions. In particular, (see ref. [20] pg. 232 (i"i7) 1;n)n+1 ), 

L 2 {a =F 2)f n ±i = (a =F 2) — {D\ t D 2 )f n±1 ■ f n = . (3.28) 

The DiS are operators of the Hirota calculus which are defined in ref. [20] for example. 
Therefore, we get 

L m (a)f n = L m (a T 2)f n±1 = , m > -1 . (3.29) 

If we now consider p = a/\^2 as the zero-mode, 'momentum' operator of the scalar field, 
with pf n = —y/2nf n , in accordance with (3.24), then its conjugate variable or 'position' 
operator is q = —it-i/\/2. This is deduced from equation (2.5) and (3.29). 



6 However, the 'topological' point described in ref. [9] does not seem to require any particular 
value of a. 
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The above result (3.29) also implies that the whole Hirota hierarchy admits a 'master' 
string equation. It is most suggestively written in terms of the full tau-function f , for 
which the f n are the projections onto eigenspaces of the zero-mode. The master string 
equation is the m = — 1 version of the the following Virasoro constraints 

L m f = 0, m > -1 , 

where the L m are the Virasoro generators of the bosonic field. So the final set of constraints 
are exactly analogous to those for the 1-arc KdV case, the difference being that there one 
has a twisted scalar field, whereas here we have an untwisted scalar field. 

The appearance of the parameter a is rather mysterious, since it was not manifest in 
the matrix model. It seems to label different sectors in the theory which are not connected 
by the flows. It is clearly desirable to have a better understanding of its origin and meaning. 



3.4 KDV Reduction of Virasoro Constraints 



We have already shown how the the ZS hierarchy can be reduced to the KdV hierar- 
chy, and how the string equations respect the reduction. On general grounds, one would 
anticipate that this would extend to all the Virasoro constraints, and we now prove this. 
The reduction involves taking ip = — e*- 1 and t 2k+ i — > Vfc. Then u = —ipip satisfies the 
KdV hierarchy with 

*i = 2- fc ^ dv , r zs = (r KdV ) 2 (3.30) 

and 

dxi 

Q t Kdv = ^xRk+i ■ (3.31) 

Notice that the second equation of (3.30) implies that under the reduction, i.e. on the 
subspace t 2k +i = V/c with ip = — e t_1 , Zzs — > ^Kdv, as it should. Under such reduction, 
H 2k = -H' 2k _ x : 

F 2k + G 2k = (F 2fc _! + G 2k -i)' + (g + /)# 2 fc-i = (g- f)H 2k -i 

= F 2k+1 + G 2k+1 = (F 2k + G 2k )' + (g + f)H 2k (3.32) 
^((g + f)H 2k . l )' + (g + f)H 2k = 

but (g + f)' = 2ip' = 0, and the result follows. 

We now show that the ZS Virasoro constraints on r zs (= f n ) reduce to Virasoro 
constraints on r KdV (= r ) for a precise value of the zero-mode a (or n). Let us first 
consider the equation corresponding to Iq, under the reduction 

J2(2k + l)t 2k H 2k + a = - J2(2k + l)* 2 fc#2fc-i + « = (3.33) 

k>l k>l 

Using the relations with the correlation functions, and integrating twice over to = we 
get 

J2(2k + l)t 2k (0 2k _ 2 ) -^i=0 (3.34) 

k>i 
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which will produce the reduced L_i constraint. The reduction of I\ is direct, and leads to 

+ l)t 2k (0 2k ) + ^ = (3.35) 

A:>0 

which produces the Lq constraint. This constraint is also obtained from I 2 , the L\ con- 
straint in the ZS hierarchy. The corresponding equation is 

J2(k + l)tkH' k+2 + aH[ + 2H 2 = 

k>0 

(3.36) 

J2( 2k + ^^fc^fc+i + ( 2 - <*) H i = 

fc>0 

Using the relations with the correlation functions, and integrating three times over to = 
we get 



J2(2k + l)t 2k (0 2k ) + (l + a)F 



k>0 



= (3.37) 



which, by consistency with eq. (3.35), requires a = —1. Therefore, the KdV reduction is 
only consistent for this value of the, a priori, arbitrary parameter a. 

Let us now consider the equation corresponding to ^3, which, again, will to fix the 
function on which the constraints act: 

^2(2k + l)t2fc#2fc+3 + MPPP) (P) + MPP) 2 + 2H 3 + 2(1- a) (PPPP) = . (3.38) 

fc>0 

In the usual way, we get 

^(2A; + l)t 2k (0 2k+2 ) + (P) 2 + (1 - a){PP) = . (3.39) 

fc>0 

This equation will lead to the L\ constraint, and, again, the second order terms fix the 
functional on which the Virasoro constraint act. In this case a = — 1, and 

4 



(P) 2 + (1 - a)(PP) = (P) 2 + 2(PP) = ^==dlJz^ ^ ZKdv^r^ . (3.40) 



KdV 



Notice that the required value of a is consistent with the quantization proposed in (3.24). 
Therefore, the reduced Virasoro constraints act on the square root of the partition function, 
which is the tau- function of the KdV hierarchy, in agreement with [11]. 

In terms of the variables t k = tf dV = 2 fc t|f , the reduced Virasoro constraints are 
L„r KdV = 0, with n > —1, where the operators L n are those of eq. (3.4). 
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3.5 MKDV Reduction 



In the case of the mKdV reduction we have already shown in section §2.2 that G 2 k+ii 
F 2k and H 2k vanish. It is straightforward to verify that under the reduction, i.e. on the 
subspace t 2k +\ = V7c with ip = e 2t_1, i/>, Zzs — > 2 m Kdv- Using these results one can apply 
the reduction directly on the constraints. 

From I (L_i constraint), we get 

^2(2k + l)t 2k H 2k + a = a = 0. (3.41) 

k>i 

Therefore, the mKdV reduction requires a = 0, which is clearly consistent with the quan- 
tization of a proposed in (3.24). From Ii, we get 

^(2£; + l)t 2k H' 2k+1 + 2H 1 = 0. (3.42) 

A;>0 

But, 

dF" 

H' 2k+1 = 2{PP0 2k ) = -2 — (3.43) 

and we get the equation 

dF" 

^(2fc + l)t2k^— + 2F" = (3.44) 

fc>0 2fc 

which has the form of an L constraint. This equation can also be rewritten as 

Y^(2k + l)t 2k -^ + v = (3.45) 

which, after integration, becomes the string equation of the mKdV hierarchy, eq. (2.20). 

If we reduce the Lq constraint of the ZS hierarchy itself, we find that the partition of 
the mKdV model satisfies 

£(2* + l)C KdV ^Kdv^mKdv = . (3.46) 

fc>0 ° l k 

Notice that fi, the parameter of the mKdV Virasoro constraints of §3.1, is determined by 
the reduction to be j^. 

In a similar way, from I 2k we get equations which are identically zero and do not give 
rise to any constraint. Instead, from I 3 (L 2 constraint) we get 

X>* + l)^ KdV ^v + 1^2 j ^Kdv = (3.47) 
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where £™ KdV = 2 k t^. We can write the above constraint in the following way. Firstly, we 
express the partition function in terms of the tau- functions 2 m Kdv = toTl- Then we use 
the relation (2.25) to write 

c^ron) = 27tfn + 2t t[', (3.48) 
from which we deduce that (3.47) may be rewritten as 

(LiTo) n + r (Lin) = , (3.49) 

which is clearly a consequence of the L\ constraints for tq and n that we found for the 
mKdV model in §3.1. 

One could carry on this process of reducing the higher Virasoro constraints. The 
resulting constraints would act directly on the partition function of the mKdV model; and 
hence would not be Virasoro constraints. Nevertheless, we expect that the constraints 
on the partition function should be expressible in terms of Virasoro constraints acting on 
each tau-function separately, as we found for above for L\. In any case, one would always 
find an expression which was compatible with the mKdV Virasoro constraints: this is due 
to the fact that the mKdV Virasoro constraints follow directly from the mKdV string 
equation and the recursion relations of the hierarchy, which are obtained by the reduction 
from the ZS string equation and hierarchy. 

4. Discussion and Open Problems 

In this paper we have attempted to analyse all the possible double scaling limits of 
the hermitian and anti-hermitian 1-matrix models. As it is clear from the fact that eq. 
(1.3), after the double scaling limit, gives rise to a differential operator of degree two in x, 
the hermitian and anti-hermitian 1-matrix models are related to the si (2, C) hierarchies. 
Hermitian and anti-hermitian matrix models have many common properties: in the 1-arc 
sector they both give rise to a KdV hierarchy for an even potential and to a doubled KdV 
hierarchy for a general potential; in the 2-arc sector, with even potential, they both give 
rise to the mKdV hierarchy. Instead they differ in the 2-arc sector with a general potential 
where the hermitian models give rise to only half of the critical points associated to the 
NLS hierarchy whereas the anti-hermitian model gives all the critical points associated to 
the ZS hierarchy (except for the 'Topological' one). 

For the hermitian models the multi-critical points obtained, although described by a 
NLS hierarchy, actually have solutions which are described by a mKdV hierarchy. We do 
not know whether these solutions are the only ones. Furthermore, these critical points 
correspond to purely even potentials; in this sense the NLS structure is irrelevant, and one 
is really dealing with a mKdV structure. Instead the ZS hierarchy admits also a reduction 
to KdV, in other words in the 2-arc sector of the anti-hermitian matrix model we found a 
new series of multi-critical points described by the KdV hierarchy, besides the one already 
known from the 1-arc sector. This set of KdV multi-critical points are not in any simple 
way connected with those in the 1-arc sector, since, for example, the topological critical 
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point describing topological gravity [14] , which cannot be obtained from the 1-arc sector, is 
obtained from the 2-arc sector with anti-hermitian matrices with a fourth order potential. 

The situation in the 2-arc sector with anti-hermitian matrices and a general potential 
seems to be the richest, being described by the ZS hierarchy. The 'even' multi-critical 
points admit solutions described by KdV and mKdV hierarchies, which require two par- 
ticular values of the parameter a. Clearly, it would desirable to understand the role of 
the parameter a, from the point of view of the matrix model, and also to know whether 
the KdV and mKdV solutions exhaust the possible solutions for the 'even' multi-critical 
points. For instance, are there other solutions for different values of a, and do solutions 
exist only for the discrete values suggested by the tau-function formalism? An interest- 
ing open question regards the nature of the 'odd' multi-critical points of the ZS string 
equation. It is now well-known that solutions for the KdV and mKdV systems describ- 
ing multi-critical behaviour exist [1,12,24]; we do not have any arguments to show that 
solutions can be found for the 'odd' scaling points of the ZS hierarchy, except for the first 
(or 'topological') point, corresponding to ti 7^ 0, which was investigated in [9]. This point 
cannot be obtained from the matrix model and, as described in ref. [9], could give rise 
to a new kind of 'topological' theory. Anyway, since the mathematical apparatus exists 
for tackling the issue of the existence of solutions to these non-linear differential equations 
[15], we hope these questions will be addressed and solved elsewhere. 

One of the results of this paper is the realization of the rather universal nature of the 
Virasoro constraints and the fact that they act on the tau-functions of the appropriate 
hierarchy, and not, necessarily, directly on the partition function. For the KdV and ZS 
hierarchies, the Virasoro constraints are 

L n r = n > -1 , (4.1) 

where in both cases r is an element of the basic representation of the Kac-Moody algebra 
for the twisted and untwisted constructions, respectively. In the mKdV case the 
partition function is the product of two tau-functions, which arise from the two basic 
representations of A± , which both satisfy the Virasoro constraints of a highest weight 
vector: 

L n Tj = fJLTj S nfi , (4.2) 

for j = and 1. 

Near the completion of this work our attention was drawn to refs. [25,26] which discuss 
similar Virasoro constraints to those above. 
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